Contractile units in disordered actomyosin bundles arise from F-actin buckling 



Martin Lenz 1 , Todd Thoresen 2 , Margaret L. Gardel 1,2,3 , and Aaron R. Dinner 1 ' 2 ' 4 
1 James Franck Institute, 2 Institute for Biophysical Dynamics, 3 Department of Physics, 
4 Department of Chemistry, University of Chicago, Chicago IL 60637, USA 

Bundles of filaments and motors are central to contractility in cells. The classic example is stri- 
ated muscle, where actomyosin contractility is mediated by highly organized sarcomeres which act 
as fundamental contractile units. However, many contractile bundles in vivo and in vitro lack sar- 
comere organization. Here we propose a model for how contractility can arise in actomyosin bundles 
without sarcomeric organization and validate its predictions with experiments on a reconstituted 
system. In the model, internal stresses in frustrated arrangements of motors with diverse veloci- 
ties cause filaments to buckle, leading to overall shortening. We describe the onset of buckling in 
the presence of stochastic actin-myosin detachment and predict that buckling-induced contraction 
occurs in an intermediate range of motor densities. We then calculate the size of the "contractile 
units" associated with this process. Consistent with these results, our reconstituted actomyosin 
bundles contract at relatively high motor density, and we observe buckling at the predicted length 
scale. 



Contractility arising from interactions between myosin 
molecular motors and actin filaments (F-actin) is used 
ubiquitously by cells to build tension and drive morpho- 
logical changes [1 . Such force transmission from molec- 
ular to cellular length scales is well understood in stri- 
ated muscle, where it critically relies on highly organized 
structures known as sarcomeres [Fig.[]Ja)] 0. However, 
many contractile actomyosin bundles found in vivo, such 
as smooth muscle fibers [3], graded polarity bundles [4] 
and the contractile ring [5], lack a sarcomeric organiza- 
tion. Most recently, we have shown that in vitro bundles 
lacking apparent sarcomeric organization can also con- 
tract [6 [e.g., Fig. [TJb)] . In these disparate systems, 
contraction occurs with a well-defined contraction veloc- 
ity per unit length, suggesting that contractile bundles 
can be meaningfully divided into elementary units that 
are arranged in series [5H?]. The mechanisms giving rise 
to such units in the absence of sarcomeric organization 
are not understood. 

Much theoretical work on non-sarcomeric actomyosin 
assemblies posits contractility as a fundamental assump- 
tion, and predicts larger-scale effects such as polarity or- 
ganization [8 , the appearance of topological defects [9], 
active stiffening [10] , and oscillatory behavior in cells [IT] . 
Models that address the microscopic origin of contractil- 
ity assume that myosin motors dwell at the barbed ends 
of F-actin, thus acting as transient static cross-linkers 
[T2] . This generates sufficient sarcomere-like organiza- 
tion to elicit contraction [13] . Experimental evidence for 
this behavior is unfortunately lacking [14] , and it is thus 
important to investigate alternative mechanisms. 

Here we demonstrate another route to contractility 
in non-sarcomeric bundles through theory and exper- 
iments. We first show experimentally that contrac- 
tion in reconstituted actomyosin bundles is accompa- 
nied by F-actin buckling. Such buckling could provide 
a symmetry-breaking mechanism necessary for contrac- 
tion [15] [Fig. flic)], and we investigate the general conse- 



quences of asymmetric filament response theoretically by 
considering the build-up of forces in a bundle with ran- 
domly arranged motors. We predict that buckling, and 
thus contraction, can only be achieved in an intermedi- 
ate range of motor density, due to the stabilizing effects 
of strong cross-linking at high density and of stochastic 
motor detachment at low density. This picture yields 
a characteristic length scale between two buckles, which 
provides a natural size for a contractile unit. These pre- 
dictions are consistent with experimental observations, 
suggesting that buckling is necessary for contractility in 
non-sarcomeric actomyosin bundles. 

To form reconstituted actomyosin bundles, we incu- 
bate F-actin with length £f — 5 jam with smooth mus- 
cle myosin thick filaments of length ~ 300 nm in buffer 
lacking ATP such that thick filaments cross-link F-actin 
with high affinity. While flexible motors have been con- 
sidered as a basis for contraction [16] , this is unlikely 
to apply here as thick filaments are significantly more 
rigid than F-actin. The bundle lengths range from 10 to 
100 jam with 4-6 F-actin per bundle cross-section, and no 
sarcomeric organization is observed [6]. By varying the 
concentration of myosin filaments, the average spacing 
£o between two consecutive myosin filaments can be var- 
ied from 390 nm to 5.3 jam [6 . Perfusing buffer contain- 
ing 1 mM ATP causes bundles formed with high myosin 
density (£q — 540 nm) to shorten by ~ 10% rapidly 
(100-600 nm-s- 1 ) [Fig.Qd) and Ref. [14 ; Movie SI]. In 
contrast, contraction does not occur at low myosin den- 
sity (£q = 1.5 jam) [Fig. [TJe) and Movie SI]. A sharp 
transition between those two behaviors is observed at 
4 = l-3/im [Fig.^f)]. 

In considering these observations, it is important to 
recognize that actomyosin interactions can a priori elicit 
extension just as well as contraction. As shown in 
Fig. [TJc), elementary bundles comprised of two F-actin 
("filaments") and one myosin thick filament ("motor") 
containing numerous myosin heads contract when the 
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FIG. 1. Contraction in actomyosin bundles, (a) Sarcomeric 
structure as in striated muscle, where passive cross-linkers 
are attached near F-actin barbed ends. Motors move towards 
F-actin barbed ends, causing each contractile unit (sarcom- 
ere) to contract, (b) Bundle devoid of sarcomeric organiza- 
tion or passive cross-linkers, as in our experiments, (c) Mo- 
tors and polar filaments induce local contraction or exten- 
sion depending on the geometry of their assembly. Overall 
contractility requires breaking the balance between these two 
effects, (d) Time-lapse images of a bundle comprised of F- 
actin and fluorescent myosin thick filaments (inverted con- 
trast) with £o — 540 nm. The initially wavy bundle becomes 
taut following the addition of 1 mM ATP at t = s, indicat- 
ing contraction. Scale bar, 5 /mi. (e) Similar experiment with 
£o = 1.5 /mi, showing no contraction. Scale bar as in (a). See 
also Movie SI. (f) Bundle contraction as a function of £o. 
Bars indicate standard deviation (n ^ 25). 



motor is located in the vicinity of the pointed ends, but 
extend when it is close to the barbed ends. In sarcom- 
eres, myosin is restricted towards the pointed ends of 
actin filaments, which favors contraction over extension 
[Fig. [TJa)]. By contrast, we show in Refs. p~4j [15] that 
these two tendencies compensate and prevent overall con- 
traction in non-sarcomeric bundles unless two specific 
conditions are fulfilled. First, the unloaded motor veloc- 
ities need to have a certain dispersion among the motor 
population. Otherwise, motors merely induce filament 
translation without inducing overall bundle contraction 
[15] . Such dispersion has been observed experimentally 
[T7] . and in our system likely arises from the variation 
of number of myosin heads in the thick filaments. In its 
presence, stresses build in the bundle as motors translate 
the filaments with different preferred velocities. While 
these stresses present a potential for bundle deforma- 
tion, the disordered nature of the bundle implies that 
they have equal chances of being compressive or extensile. 
Thus the second condition is an asymmetric response of 
the filaments to such stresses, which breaks the symmetry 
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FIG. 2. Buckling in non-sarcomeric contractile actomyosin 
bundles, (a) Time-lapse images of fluorescent actin (inverted 
contrast) showing F-actin buckling (arrowheads) following the 
addition of 1 mM ATP at t = s. Scale bar, 5 /im. See 
also Movie S2. (b) Relative contraction (filled squares) and 
number of F-actin buckles (open circles) as a function of 
time. Data shows mean ± sd averaged over n — 3 bun- 
dles with £o c± 1 /mi. (c) The presence of fast (grey) and 
slow (white) motors generically induce compressive (red) and 
extensile (blue) stresses in filaments, (d) Buckling of the com- 
pressed filaments leads to an overall shortening of the bundle. 



between contraction and extension. Contraction occurs 
when filaments yield under compression while resisting 
extension [14] [15] . 

We next look for evidence of this behavior in our ex- 
periments. In this respect, we observe F-actin buckling 
coinciding with contraction [Fig. [2ja); Movie S2]. This 
constitutes an extreme form of asymmetric response of 
the filaments, and thus enables contractility. Prior to 
ATP addition, compact bundles with aligned F-actin are 
observed. Upon ATP addition, the frequency of buckles 
increases rapidly during contraction, and then diminishes 
once contraction stops [Fig. |2jb)]. These F-actin buckles 
are dynamic, with their amplitude, curvature and loca- 
tion changing over time. 

Qualitatively, the relationship between buckling and 
contraction can be understood as follows. Consider two 
antiparallel filaments interacting through several differ- 
ent motors with distinct speeds [Fig. |2jc)]. As motors 
start to move relative to the filaments, stresses build in 
sections of the filament flanked by motors with differ- 
ent speeds. When the flanking motor proximal to the 
barbed end is faster than that proximal to the pointed 
end, compression arises. When it is slower, tension arises. 
Following buckling of the compressed filament sections, 
fast motors are free to move quickly while the others move 
slowly. This results in the growth of the compressed sec- 
tions and shrinkage of the extended ones, and thus in 
overall bundle contraction [Fig. [2jd)]. The region cen- 
tered around each buckle thus plays the role of a con- 
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tractile unit, whose typical size is equal to the distance 
£b between two buckles. 

In this picture, the contractile behavior of the bundle 
hinges on the ability of the motors to induce filament 
buckling. At high motor density, we expect the bundle 
to be so strongly cross-linked that buckling becomes im- 
possible despite the sizable stresses induced by a large 
number of motors. At low motor density, we expect that 
stochastic detachment of the motors undermines stress 
build-up and thus prevents buckling. Here we present a 
mathematical model to predict the range of myosin den- 
sities enabling contraction and the contractile unit length 
£b- These results are then compared with the observa- 
tions in Figs, [l] and [2] to validate the proposed contraction 
mechanism. 

The key assumptions of our model are that (1) mo- 
tors have a dispersion in their unloaded velocities, (2) a 
section of filament between two motors buckles above a 
certain threshold force Fb, and (3) motors intermittently 
detach from the filaments, thus allowing local stress re- 
laxation. We consider a bundle of weakly deformed fil- 
aments and ask whether the forces developing within it 
are sufficient to induce buckling [Fig. [3ja)]. 

To this end, we focus on a single filament of length if 
and approximate its surroundings by an effective medium 
composed of immobile [14 point-like motors separated by 
a distance 1$ <C if [Fig. [3^b)]. This divides the filament 
into discrete sections, which we label by i = 0, . . . , if/io- 
We take into account the possibility that the filaments 
are not straight, but bend away from the x-axis, imply- 
ing that the contour length Li of filament section i can 
be larger than i$. Defining fi as the tension of filament 
section i (fi < for a compressed filament section), we 
expand its force-extension relationship for small defor- 
mations: 



Li = Li(fi = 0) - eft, 



(1) 



where c > is the filament compliance. We refer to 
the motor flanked by filament sections i — 1 and i as 
"motor z", and describe its operation by the simplified 
force-velocity relationship 



fi-i - fi 



(2) 



Here V{ denotes the local velocity of the filament at the 
location of motor i and \ > is the motor susceptibil- 
ity. Eqs. ([T]) and Q yield a local relaxation time scale 
r r = x c /2- The time-independent stall force of motor i 
is denoted by Fi in Eq. Q , and is drawn from a random 
distribution satisfying 



Fi=F s and F % F 3 - Fi F 3 = 5F 2 s 5 l3 , 



(3) 



where bars denote averages over the motor distribution. 
As a result, different motors have different unloaded ve- 
locities Fi/x as required for contraction. Owing to the 




7777777777777777777777777777777777777777777777777777777777777 

1 2 ••• i-1 i i + 1 ••■ £ f /£ 


FIG. 3. Stress build-up in bundles with non-identical motors, 
(a) In a bundle with motors having non-identical velocities 
(shades of grey), filaments of lengths w if are subjected to 
random motor forces at points ~ £0 apart distributed through- 
out their length, (b) Prior to buckling, the environment of a 
filament of interest (red) can be approximated by a collection 
of immobile motors (shades of grey) [14] . 



conservation of filament mass: 



H+l- 



(4) 



Finally, a motor bound to several filaments as in Fig. [3ja) 
can transiently detach from one while still holding onto 
the others [14] . We thus let each motor i randomly detach 
from the filament with a constant rate 1/r^. Following 
detachment, local filament stresses relax instantaneously, 
yielding fa = fi_ x = (fi + /»_i)/2. The motor then 
reattaches after a time much shorter than r r and [14] . 
We denote by (...) the average over the Poisson process 
of motor detachment. 

We obtain the space and time evolution of the filament 
tension f(x,i) in the continuum limit i —> x/io by com- 
bining Eqs. ([l]|4| and averaging over motor detachment 

Eh 



d t (f) - Dd 2 x (f) = (£o/2r r )d x F, 



(5) 



where D = il(r~ x + r~ 1 )/2. The right-hand-side of 
Eq. ([5| involves the spatial gradient of the stall force 
F(x), reflecting the fact that non-identical motors lead 
to force build-up. This effect competes with the relax- 
ation of filament forces through motor detachment, which 
enters through the diffusion term Dd^(f). 

An initially relaxed filament [f(x, t = 0) = 0] experi- 
ences a vanishing average force {f)(x,t) = throughout 
its dynamics. To quantify the magnitude of the motor- 
induced stress, we thus calculate the rms filament force 
((Z 2 )) 1 / 2 [14 and find that it increases monotonically 
from zero at t = to 
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at t = 00 [Fig. [4ja-b)]. We next estimate the depen- 
dence of the ratio r r /rd on the experimentally accessible 
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FIG. 4. Model predictions for filament force build-up. 
(a) Black line: Steady-state filament force foo as a function 
of motor spacing £q as in Eq. (J6|. For £q <C and ^> ^o, 
/oo oc ^^ 1 ^ 2 and ^~ 9 ^ 2 , respectively. Colored lines: buckling 
force Fb oc (b) Typical filament force ({f 2 )) 1 ^ 2 as a func- 
tion of time (c) Contractile unit size £b as a function of 
as in Eq. (|8| (£/ c± 5 /im). 



merits. 

To characterize the contractile units resulting from this 
mechanism when £q < £q < £q , we turn to the transient 
regime leading up to filament buckling. The filament 
force profile as a function of x is initially flat, and sub- 
sequently coarsens into a random walk for t = +oo. Ac- 
cording to Eq. ([5| , this coarsening occurs diffusively with 
diffusion coefficient D. The typical filament forces at 
time t < £ 2 /D are thus of order foo(VDt/£f) 1/2 . We de- 
note the time that this force reaches the buckling thresh- 
old Fb by ts , following which contraction proceeds as in 
Fig. [2|c-d) and the coarsening dynamics is interrupted. 
The distance between buckles at £b thus yields the con- 
tractile unit size 



(8) 



parameter £q. A Worm-Like Chain model for filament 
elasticity yields c « £\jkBT£ 2 ^ where £ p is the filament 
persistence length [18] , and we approximate x ~ Fs/ V , 
where v is a characteristic motor velocity. This implies 
T r /r d « {£ Q /£l)\ with £* = (kBT^vrd/Fs) 1 ^. We can 
thus distinguish two regimes for the steady-state force 
foo [Fig. Qa)]. For £q <C detachment events are rare 
compared to the time r r needed for the force to recover 
from such an event, and f^ is not affected by them. For 
A) ^> Zoo quickly decreases with increasing £q as de- 
tachment becomes much faster than recovery. 

Up to a prefactor of order one, contraction proceeds 
as in Fig. d^c-d) if > F B « k B T£ p /£l [18]. Com- 
paring /oo to Fb as in Fig. |4ja), we find a threshold 
stiffness above which buckling cannot occur (as exam- 
pled by the blue line). Reasonable values for our ac- 
tomyosin system are £ p ~ 10 /im, v ~ 200 nm • s _1 , 
« ~ 1 pN and ~ 200 ms based on the typical 
time scales involved in the myosin mechanochemical cycle 
[T4] . These values put us in the soft filament regime de- 
fined by £ p < 5FlL 2 f (vT d ) 3 / 2 /k B T 5 / 2 Fs /2 ~ 20 cm (red 
line). In this regime, the lines representing Fb and f^ 
intersect at 

q = (k B T£ p /5Fs£) /2 ) 2/s ~ 70 nm, (7a) 
£+ = (£) /2 VT d £ P SF s /F s ) 2/b ~ 1 /im, (7b) 

meaning that buckling and contraction occur for £q < 
< £q • This range reflects the fact that strong cross- 
linking (£0 < £0 ) suppresses buckling while sparse mo- 
tors (£q > £q) are undermined by stochastic detachment. 
While the regime £0 ~ £q is not accessible experimen- 
tally, the predicted value for £q is strikingly similar to 
the motor spacing at which the breakdown of contrac- 
tion is observed in Fig.[TJf) (1.3 /im), suggesting that the 
proposed mechanism is a good description of our experi- 



As illustrated in Fig.[4jc), £b is typically in the microme- 
ter range, in agreement with the observations of Fig.[2|a) 
and the findings of Ref . [6] . 

Because of compensating effects between contractile 
and extensile motor-filament configurations, the familiar 
framework involving rigid filaments and identical motors 
commonly used to describe striated muscle contraction is 
not suited to study actomyosin bundles lacking sarcom- 
eric organization. Here, we put forward an alternative 
mechanism based on our observation of buckling. The 
buckling arises from the nonlinear elastic response of F- 
actin [19] and dispersion in the speeds of myosin motors 
[20] . F-actin buckling has previously been invoked to ex- 
plain contraction qualitatively [21 . Addition of passive 
cross-linkers, which are formally equivalent to immobile 
motors, would reinforce a dispersion of motor velocities 
and promote contraction. 

The order-of-magnitude agreement between theory 
and experiments with respect to the size of contractile 
units and the critical myosin concentration required for 
contraction suggests that our current analysis offers a 
good description of the onset of bundle contractility. Our 
conclusions are robust to inclusion of features such as in- 
homogeneous motor spacings £0 and force dependence of 
the motor detachment rate [14]. Our mechanism is a 
general one and applies to any one-dimensional system 
of polar filaments and motors. Further experiments and 
theory are needed to better understand the molecular ba- 
sis for motor inhomogeneities and filament asymmetric 
response in the myriad of non-sarcomeric organizations 
found in vivo. 
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